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Abstract. The Bose Hubbard dimer Hamiltonian is a simple yet effective model for descri- 
bing tunneling phenomena of Bose Einstein condensates. One of the significant mathema- 
tical properties of the model is that it can be exactly solved by Bethe ansatz methods. Here 
we review the known exact solutions, highlighting the contributions of V.B. Kuznetsov to 
this field. Two of the exact solutions arise in the context of the Quantum Inverse Scattering 
Method, while the third solution uses a differential operator realisation of the su(2) Lie 
algebra. 
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1 Introduction 

The experimental realisation of Bose-Einstein condensation using atomic alkali gases has pro- 
vided the means to study macroscopic tunneling in systems with tunable interation parame- 
ters [9]. From the theoretical perspective, the Bose-Hubbard dimer model (see equation (pQ) 
below), also known as the discrete self-trapping dimer [HO [6] or the canonical Josephson Hamil- 
tonian [9] , has been extremely useful in understanding this tunneling phenomena in the context 
of a bosonic Josephson junction. Despite its apparent simplicity, the Hamiltonian captures the 
essence of competing linear and non-linear interactions which lead to non-trivial dynamical be- 
haviour and ground-state properties (e.g. [TJ [121 Q2J EE])- I R particular the model predicts 
macroscopic self-trapping and the collapse and revival of Rabi oscillations, features which have 
been directly observed experimentally in a single bosonic Josephson junction [H 111] . 
The Bose-Hubard dimer Hamiltonian is given by 

H = k -(N l -N 2 f-^(N l -N 2 )- £ -(b\b 2 + b\b l ), (1) 

where b\, b\ denote the single-particle creation operators for two bosonic modes and N\ = b\b±, 
N 2 = b\p2 are the corresponding boson number operators. The coupling k provides the strength 
of the scattering interaction between bosons, f/, is the external potential and £ is the coupling 
for the tunneling. The change £ — > —£ corresponds to the unitary transformation b\ — > b±, 
b 2 — ► — £>2> while fx — > —[i corresponds to 61 <-> b 2 . The total boson number N = N\ + N2 is 
conserved and consequently the model is integrable as it has only two degrees of freedom and 
two conserved operators, viz. H and N. Mathematically the Hamiltonian is of interest because, 
related to its integrability, it admits exact Bethe ansatz solutions. This property opens avenues 
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to rigorously analyse the model. For example, the Bethe ansatz solution can be used to study 
the ground-state crossover from a delocalised state to a "Schrodinger cat" state in the attractive 
case [3] , as well as facilitating the calculation of form factors [10] . 

The first Bethe ansatz solution of the Hamiltonian was given by Enol'skii et al. [H [5] using 
the machinery of the Quantum Inverse Scattering Method. A key ingredient in this approach 
was the use of a bosonic realisation of the Yang-Baxter algbera, which was developed in the 
work of Kuznetsov and Tsiganov [8]. For zero external potential an alternative application of 
the Quantum Inverse Scattering Method, using the Gaudin algebra formulation, was given by 
Enol'skii, Kuznetsov and Salerno [6]. We remark this method of solution for the model has also 
been recently discussed in |13[ [Li] . In this approach a connection was made with confluent Heun 
polynomials. It was also observed in their work [6] that this connection could be established 
using an su(2) realisation of the Hamiltonian (see also p2]). This property provides a direct 
route to a third Bethe ansatz solution using elementary properties of second-order ordinary 
differential eigenvalue equations with polynomial solutions. 



2 Exact Bethe ansatz solution I 



In this section we review the Quantum Inverse Scattering Method and associated algebraic Bethe 
ansatz. The notational conventions we adopt follow those of [10], which also contains the full 
details for the following calculations. Then we will apply this approach to derive the exact Bethe 
ansatz solution of ([1]), as was originally described in [HE]. 

We begin with the su(2)-invariant i?-matrix R(u) £ End(C 2 <S>C 2 ), depending on the spectral 
parameter mEC: 



R(u) 
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(2) 



with b(u) = u/(u + rf) and c(u) = r]/(u + r]). Above, r] is an arbitrary parameter. It is easy to 
check that R(u) satisfies the Yang-Baxter equation 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v) (3) 

on End(C 2 ®C 2 ®C 2 ). Above Rjk(u) denotes the matrix acting non-trivially on the j-th and 
k-th spaces and as the identity on the remaining space. Next we define the Yang-Baxter algebra 
with monodromy matrix T(u), 



T(u) 



A{u) B{u) 
C{u) D{u) 



subject to the constraint 

Ri 2 (u - v)Ti(u)T 2 (v) = T 2 (v)Ti(u)Ri 2 (u - v). 
Given a representation ir of the monodromy matrix, the transfer matrix is defined 

t(u) = ir(A(u) + D(u)) 



(4) 



(5) 



(6) 



which satisfies [t(u),t(v)] = for any choice of u and u as a result of ([3]). If there exists 
a pseudovacuum state \x) which satisfies 



7T(A(u))\x) = a(u)\ X ), n(B(u))\ X )=0, 
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tt(C(u)) |x> ^ 0, 7r(Z)(u))| X > = d(«)|x) 
the transfer matrix has eigenvalues 

M M 
A(u)=o(u)M +d(u)\\ . (7) 

- LJ - U — Vk - LJ - U — Vk 

k=l fc=l 
Provided the Bethe ansatz equations 

^4=n^^, *=i m (8) 

are satisfied. 

We may choose the following realization for the Yang-Baxter algebra, with arbitrary weC, 

n{T(u))=L\{u + u)L h 2 (u-u) (9) 

written in terms of the bosonic realisation of the Lax operator given by Kuznetsov and Tsiga- 
nov [8]: 

u + rjNi 

V 



L»=[ 6 t' Ji), * = 1,2. (10) 



Since L(it) satisfies the relation 

12i2(u-«)Lj 1 («)L^(t;) = Lj 2 (t;)Lj l («)12 12 («-t;), i = l,2 (11) 

it is easy to check that the relations of the Yang-Baxter algebra ([5]) are obeyed. Specifically, 
the realisation of the generators of the Yang-Baxter algebra is 

n(A(u)) = (u 2 - lo 2 )I + T]uN + rj 2 N 1 N 2 - ^(Nx - N 2 ) + ftjbi, 
tt(B(u)) = (u + uj + r]Ni)b 2 + rf x b x , 
ir(C(u)) = b\{u - u + V N 2 ) + rj^bl 
n(D(u)) = b\b 2 + tT 2 /, 

It is straightforward to verify the Hamiltonian ([1]) is related with the transfer matrix ([6]) through 

1 

where the following identification has been made for the coupling constants: 

k prf n £ 
8=—' 2=- p ^ 2 =P - 

We can apply the algebraic Bethe ansatz method, using the Fock vacuum |0) as the pseudova- 
cuum \x), giving 

a(u) = u 2 — u 2 , d(u) = r]~ 2 . 

For this case the Bethe ansatz equations are 

M 



H = -p[ t{u) - 7(t'(0)) 2 - ut'iO) - rf 2 + uj 2 -u 2 
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where M is the eigenvalue of the total number operator N. The energies of the Hamiltonian are 

V 2 M 2 



4 

M 



UT]M — u z 



/ \ 

fiV Vi-uJ 



i=l 

This last expression is independent of the spectral parameter u, which can be chosen arbitrarily. 

3 Exact Bethe ansatz solution II 

The second Bethe ansatz solution of ((TJ described by Enol'skii, Kuznetsov and Salerno [B] applies 
only when fJ, = 0, i.e. for the Hamiltonian 

H = ^{N l -N 2 f- £ -(b\b 2 + b\b l ). (13) 
To obtain this solution, first we introduce new operators through a transformation 

h = -^( a i - * a 2)> b \ = ( a i + «4)» 

b 2 = -^{ai+ia 2 ), b\ = -^=(a\ - ia\) 

such that the canonical commutation relations [aj,ajj = Sj^I etc. hold. Under the above trans- 
formation the Hamiltonian (1131) becomes 



H = - I -(2m + /)(2n 2 + I) - (aj)^ - " ^ ) + 3 (n2 " m) ' (14) 

where rij = ajaj and N = n\ + n 2 . 

The next step is to write (j!4j) in terms of an sit(2) realisation. The su(2) algebra has 
generators {<S' 2 ,S' : ' = } with relations 

[5 Z ,5 ± ] =±S ± , [5 + ,5"] = 25 2 . (15) 

It may be shown that 

5+ = -i(at) 2 , ^ = ^« 2 , 5 z = i(2iV + I) 
is an su(2) realisation preserving the commutation relations (|15[) . It follows that we may write 

H = ~ + 5+5 2 " + Sr5+ - i A + S (SI - Sf) . (16) 

To derive the Bethe ansatz solution for (1161) . one takes 



exp(?/a) 
I exp(— rja) 

with a £ C, and constructs the monodromy matrix 

7T(T(u))=gLf(u + P)Ll(u-P), 
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where (5 € C and 



Lt[u >~ u V r)St ul-rjS? )> l ~ 1 ' 2 - 
The elements of the monodromy matrix are found to be 

MAiu)) = exp(„ Q ) { (/ + -l^flf) (/ + + ^S^S}} , 

»(*<„)) = e*p(„ a ) j-^Sf (/ - ^Sj) + ^ (/ + ^Sf) } , 

»(«(„)) = exp(-, a) { (/ - -S-3l) (l - ^Sj) + ^Sfflf} 

from which we can construct the transfer matrix ([6]). For the Bethe ansatz solution, the pseudo- 
vacuum state |x) can be chosen to be the vacuum state |0), either of the one-particle states 
a\ |0) or a 2 |0), or the two particle state a\a\ |0), since for all cases 

n(B(u)) \ x ) = 

and 

n(A(u)) \ X ) = a(u) \ X ) , n(D(u)) \ x ) = d{u) \ x ) . 

In this manner the form of the Bethe ansatz solution depends on whether the total particle 
number is even or odd. We find 

oW= exp(, a ,( 1 + -^)( 1 + -^), (17) 
d(«) = exp(-„ a ) ( 1 - -^p-] ( 1 - ■ (18) 



u + (5 J \ u — P 

where K\ = K2 = 1/4 or k\ = k>2 = 3/4 for the even case, and K\ = 3/4, k>2 = 1/4 or K\ = 1/4, 
K2 = 3/4 for the odd case. It can now be shown that t\, T2 defined by 

n = lim lim (^P) t(u) = 2aSf - ± (2SfS$ + S+S£ + S£S+) , 
t 2 = lim lim (^f) t(u) = 2aSI + ±- {2SfS z 2 + S+S^ + S^S+) 

are related to the Hamiltonian (|16j) and the total number operator through 
k 2 

h = T 2 -n- —I, n = -{t 1 + t 2 )-i 

lb t 

with 

To make the Bethe ansatz solution of the model explicit it is a matter of substituting (fTT|) , (fl~8|) 
into ([8]) and taking the limit ij — > to obtain 



M 

«1 K2 

a H t; + 



^ + /3 - /3 ^rivj -v k 

3r K 



M 

, k = l,...,M. (19) 
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Letting \j denote the eigenvalue of Tj , it follows from © that 

/ u + P\ I K 2 1 | 

hm lim — = — A(u = 2ki la — > , 

lim Mm ( ^--i) A(u) = 2K 2 \a + ^-- V ) . 

v-+Q^p\ v 2 J VJ I 2/3 j^Vj-P] 

The eigenvalues of the Hamiltonian are given by 

£ = £ (K2 - «l) + — + -y 2^ U i' 

J=l 

where those of the number operator are 
N = 2M + 2(«i + K2) - 1. 

It is apparent that the Bethe ansatz equations (|12|) with a; = 0, which are in multiplicative 
form, take on a different form to those given by (|19p which are additive. Moreover, the Bethe 
ansatz equations ([T2]) are associated with a single reference state whereas (fH?|) are dependent on 
the choice of reference state. In this latter case there are four forms of the Bethe ansatz equations 
associated with the choices of k±, k 2 which can take values 1/4 or 3/4. In the following it will 
be shown how a unified system of Bethe ansatz equations can be derived in the additive form. 
This approach does not use the Quantum Inverse Scattering Method. 



4 Exact Bethe ansatz solution III 

We again follow the work of Enol'skii, Kuznetsov and Salerno [6] (see also [17] ) and start with 
the Jordan-Schwinger realisation of the su{2) algebra (115|> : 

S+ = b[h, S- = blh, S' 2 = ^(iVi-iV 2 ) 

which is (N + l)-dimensional when the constraint of fixed particle number N = N± + N 2 is 
imposed. In terms of this realisation the Hamiltonian (JTJ) may be written as 

H = ^{S Z ?-^S Z -^(S + + S-). (20) 

The same (N + l)-dimensional representation of su(2) is given by the mapping to differential 
operators 

oz d N c+ AT 2 d o- d 

S =u- — , S + = Nu-u — , S = — 

au 2 au au 

acting on the (N + l)-dimensional space of polynomials with basis {l,ii, u 2 , . . . ,u N }. We can 
then equivalently represent (]20p as the second-order differential operator 



rr k ( 2 d 2 Art d N 2 \ ( d N\ £ / ^ 2x d 

ku 2 d 2 1 . . , „ . , NN d kN 2 aN £Nu 
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Figure 1. Energy levels E versus coupling k of the Hamiltonian JI])for7V=10,/i = 0, and £ = 1. 

Solving for the spectrum of the Hamiltonian (JTJ) is then equivalent to solving the eigenvalue 
equation 



HQ(u) = EQ{u) 



(22) 



where H is given by (12 If) and Q(u) is a polynomial function of u of order N. 

From this point, it is little effort to obtain a third Bethe ansatz solution for the Hamilto- 
nian (pQ) (cf. [3]). First express Q(u) in terms of its roots {vj}: 



N 



Q( u ) = Y[( u - v j)- 

Evaluating ([22]) at u = vi for each I leads to the set of Bethe ansatz equations 



Svf + (k(l - N) - 2fi)vi - £ 



kvf 



N 

E 



Vj - VI 



1 = 1, 



,N. 



(23) 



Writing the asymptotic expansion Q(u) ~ u N — u N 1 ^ Vj and by considering the terms of 
order N in (|22p . the energy eigenvalues are found to be 



N 



„ kN 2 uN ^ 

E = — + - > ' 

8 2 2 ^ 



(24) 



In the above manner a single form of additive Bethe ansatz equations (|23p is obtained. As far as 
we are aware, the mapping of the solution (J23]), (|24p to (|12p . ([2]) remains an unsolved problem. 

Fig. Q] shows the energy levels for the model with fi = and N = 10, obtained from the 
solution (|23p . (|24p . Even for such low particle number it is clearly seen the ground state be- 
comes quasi- degenerate in the attractive regime. This property underlies the validity of using 
spontaneous symmetry breaking based on a mean-field approximation, as discussed in [2], to 
distinguish the quantum phases of the Bose-Hubbard dimer Hamiltonian ([1]). Alternatively, 
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associated with the Bethe ansatz solution (|23p . (|24p there is a mapping of the spectrum of the 
Hamiltonian into the low energy spectrum of a one-dimensional Schrodinger equation. This 
facilitates a different approach for determining quantum phases of the Hamiltonian where the 
crossover is identified with a bifurcation of the Schrodinger equation potential [3]. 
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